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ABSTRACT 



Fiber optic communications are the future of U.S. Navy shipboard 
communications. They present tremendous bandwidths with no susecptability to electro- 
magnetic interference (EMI) and outstanding signal-to-noise ratios. Current technology 
uses wavelength division multiplexing (WDM) to allow multiple users on a single fiber 
simultaneously. The optical filters necessary to demodulate the WDM signal can be 
expensive. A less costly alternative could be achieved if Hadamard-Walsh Code Shift 
Keyed (CSK) encoded signals were used. An optical receiver was proposed, using the 
Mach-Zender coupler, capable of demodulating a Hadamard-Walsh CSK encoded signal. 
Building on previous work studying the characteristics of the Mach-Zender coupler, a 
relationship between the probablility of bit error and signal-to-noise ratio (SNR) was 
developed for a single user and a DPSK optical signal. This relationship was then used to 
develop an understanding of the bit-error rate to SNR relationship for a multiple-user 
CSK optical signal. Using the theoretical performance as a guide, a MATLAB model was 
then constructed to investigate the sensitivity of the receiver to non-ideal components. 



v 




VI 




TABLE OF CONTENTS 



I. INTRODUCTION 1 

H. DEVELOPMENT OF A MATHEMATICAL MODEL FOR THE MACH- 

ZENDER COUPLER 5 

A. OVERVIEW OF MACH-ZENDER COUPLER 5 

B. DEVELOPING A MATHEMATICAL MODEL 6 

1. Coupler Analysis 6 

2. Delay Element Analysis 8 

3. Phase Shift Element Analysis 9 

C. INPUT SIGNAL ANALYSIS 9 

D. SUMMARY 11 

HI. DPSK DEMODULATION ANALYSIS 13 

A. INPUT DPSK SIGNAL 1 3 

B. RECEIVER ANALYSIS 14 

1 . Output of Mach-Zender Coupler 14 

2. Photodetector Analysis 14 

3. Amplifier Analysis 15 

4. Integrator Analysis 1 5 

5. Decision Variable Analysis 17 

C. PROBABILITY OF BIT ERROR 19 

D. SUMMARY 19 

IV. CSK DEMODULATION ANALYSIS 21 

A. INTRODUCTION 21 

B. INPUT CSK HADAMARD-WALSH ENCODED SIGNAL 22 

C. RECEIVER ANALYSIS 24 

1 . Output of First Stage Mach-Zender Coupler 24 

2. Output of Second Stage Mach-Zender Couplers 25 

3. Photodetector Analysis 25 

4. Amplifier Analysis 26 

5. Integrator Analysis 27 

6. Decision Variable Analysis 28 

D. PROBABILITY OF BIT ERROR 29 

E. SUMMARY 30 

V. SIMULATION AND SENSITIVITY ANALYSIS 33 

vii 



A SIMULATION 33 

1. Design Approach 33 

2. Mach-Zender Coupler Simulation Construction 33 

3. Electrical Receiver Element Simulation 34 

4. Statistical Analysis 34 

B. SENSITIVITY ANALYSIS 34 

1 . Delay Error 34 

2. Phase Shift Error 36 

C. SUMMARY ' 37 

VI. CONCLUSION 39 

A. DISCUSSION OF RESULTS 39 

B. FUTURE RESEARCH 40 

APPENDIX A. ANALYSIS OF THE OUTPUT OF TWO MACH ZENDER 

COUPLERS FOR CSK ENCODED SIGNAL 41 

APPENDIX B. MATLAB SIMULATION PROGRAMS 45 

LIST OF REFERENCES 57 

INITIAL DISTRIBUTION LIST 59 



viii 



I. INTRODUCTION 



The majority of data transported through modem U.S. Naval platforms is carried 
on copper wire or coaxial cable. As the amount of data being transferred dramatically in- 
creases, the need for a more efficient methods of transport arises. Given the tremendous 
increase in bandwidth of fiber optics, its use to transport data in naval platforms is be- 
coming standard. Not only does fiber optic communication offer a tremendous band- 
width, it offers many other desirable characteristics that include resistance to electro- 
magnetic interference (EMI) and an excellent signal-to-noise ratio that make its use on 
naval platforms extremely attractive. 

To design the information flow in naval vessels of the future, it is important to 
study methods of optical data transmission and explore ways to improve upon and exploit 
the tremendous bandwidth available. In this thesis an alternative to Wave Division Multi- 
plexing (WDM) will be developed and analyzed. This alternative relies upon the use of a 
specialized coupler, the Mach-Zender coupler, to demodulate the multiplexed optical sig- 
nal into individual signal components. 

The Mach-Zender coupler is composed of two 2x2 optical couplers, a delay line 
and a phase shift element. The signal enters the first coupler input port (port 1) as shown 
in Figure 1.1 and is equally divided between the output ports A and B. Signal a, is then 
delayed one bit period (T) and becomes signal a \ . Both signals then enter the second 
coupler are combined and then equally divided to the output ports A and B. The signal 
b 2 is then phase shifted - n 12 radians giving the outputs of the Mach-Zender coupler of 
a 2 and b 2 . 

Utilizing the outputs of the Mach-Zender coupler, a receiver can be constructed to 
demodulate a single optical Differential Phase Shift Keyed (DPSK) signal as shown in 
Figure 1.2. The outputs of the Mach-Zender coupler are followed by a photodetector 
(PD) to convert the optical signal to an electrical current by an amount proportional to the 
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s^> 




Figure 1.1 The Mach-Zender Coupler 



responsivity ( 91 ) of the photodiode. This current signal is further amplified by a transim- 
pedence amplifier of gain Z# and input to an integrator over one period (7). The energy 
of the two channels is then compared and a symbol (“1” or “0”) is chosen based on the 
channel with the greatest energy. 




Figure 1.2 Demodulation of DPSK Signal Using Mach-Zender Coupler 

By combining multiple couplers with varying delays (e.g. T and 27), a Code Shift Keyed 
(CSK) signal employing Hadamard-Walsh coding could be demodulated. This orthogo- 
nally coded signal would allow multiple stations to use a single fiber optic line simulta- 
neously as is currently being done using Wavelength Division Multiplexing (WDM). 
With WDM, however, there is a significant cost in creating very sensitive filters neces- 
sary to demodulate the signal [Ref. 1], Figure 1.3 shows a system using multiple Mach- 
Zender couplers that is capable of decoding four Hadamard-Walsh coded signals. 
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By increasing the number of Mach-Zender coupler stages, this receiver can be con- 
structed to demultiplex increasingly complex Hadamard-Walsh coded signals. 

In Chapter II, the theory behind the Mach-Zender coupler will be discussed, and 
from this theory, a mathematical model of the coupler will be developed. This model will 
include both time and frequency domains. Chapter III will use the Mach-Zender mathe- 
matical model of Chapter II to develop a single user DPSK encoded optical receiver. 
From this DPSK receiver a probability of bit-error will be determined based on the sig- 
nal-to-noise ratio at the receiver. Chapter IV will then propose a multi-user Code Shift 
Keyed (CSK) optical receiver and develop the bit error rate for the multi-user receiver. 
Chapter V will then compare the theoretical results of the bit error rate versus signal-to- 
noise ratio in the single user receiver with those obtained by a computer model of the re- 
ceiver. Chapter V will continue with the computer model to conduct a sensitivity analysis 
of the single user receiver. Finally, Chapter VI will discuss the conclusion of this thesis 
and present some areas for future research with these optical receivers. 
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II. DEVELOPMENT OF A MATHEMATICAL MODEL FOR THE MACH- 

ZENDER COUPLER 

A. OVERVIEW OF MACH-ZENDER COUPLER 

In the following analysis, a mathematical description of the Mach-Zender coupler, 
as has been previously developed at the Naval Postgraduate School [Ref. 2], will be re- 
viewed. This model will be used to analyze the demodulation of an optical DPSK signal 
using a single Mach-Zender coupler. 

As previously described, the Mach-Zender Coupler consists of two 2x2 optical 
couplers, a delay line of some delay nT where n is an integer and T is the period of the 
signal, and a device to provide a phase shift of -n / 2 radians. Port 1 of the first coupler 
is the input and Port A of the second coupler and the output port of the final phase shift 
element are the outputs. 

The delay of the Mach-Zender coupler can be accomplished by inserting a length 
of fiber equivalent to the desired delay. Heinbaugh [Ref. 2] proposed using a 2x2 coupler 
with a single input and a single output to obtain a 7t/2 phase shift. This mathematical 
model will use a - 7 t /2 phase shift. 

The input signal is defined as s(t) in the time domain and S(f) is its Fourier trans- 
form in the frequency domain. As the review in this chapter will show, the output in the 
time domain is s(t-T)-s(t) for port A of the second coupler and s(t-T)+s(t) for the output 
of the phase shift element. This is as illustrated in Figure 2.1. 




Figure 2.1 Signal Throughput of the Mach-Zender Coupler 
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B. DEVELOPING A MATHEMATICAL MODEL 



In this section each element of the Mach-Zender coupler will be analyzed mathe- 
matically and an expression relating to that particular element’s effect on s(t) will be de- 
veloped. 

1. Coupler Analysis 

The optical coupler, Figure 2.2, has four ports labeled 1, 2, A and B. Each port 
can be used as an input or an output . 



a l~» 


1 A 


<-a 3 


<- t>! 


b 3 -> 




2x2 




a 2 -> 


2 B 


<- a 4 


b 2 


b 4 -> 



Figure 2.2 Optical 2x2 Coupler 

Let the frequency domain inputs of the 2x2 coupler of Figure 2.2 be described by 
the following vectors a and b : 




( 2 . 1 ) 



The input vector a and the output vector b are related by the scattering matrix S 
defined as follows [Ref. 3]: 
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Equation 2.2 can be simplified if the matrices of Equation 2.4 are partitioned as 
shown below: 
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After partitioning. Equation 2.5 can be rewritten as follows: 
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Observe that in the ideal 2x2 optical coupler, no part of the input signal is re- 
flected back to the input ports. Specifically, when the only inputs are a, and a 2 , the out- 
puts are confined to b 2 and as in Equation 2.6b. That is b 2 and b 4 consist of only in- 
phase and 7 r /2 phase-shifted components of o, and a 2 . 




Figure 2.3 Single Input 2x2 Coupler 

To apply this result to a practical use of the Mach-Zender coupler, consider Figure 
2.3 in which one input signal is applied to one input port, and the remaining port is termi- 
nated. Equation 2.6b now reduces to become: 
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We see that the input S', is decreased by a factor of l/ a/ 2 at ports A and B and the output 
at port B is shifted in phase by n / 2 radians from the input. 

2. Delay Element Analysis 

As shown in Figure 1.1 the delay element represents a time shift of T in the time 
domain which corresponds to multiplication by e jlnfr in the frequency domain. The input 
and output of the delay element can be modeled in the time domain as in Figure 2.4 with 
the output a I = S',e /2 ' l/7 in the frequency domain. 
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Figure 2.4 Time delay element 



3. Phase Shift Element Analysis 

The phase shift of the Mach-Zender Coupler can be accomplished using a 2x2 
coupler with only one input and one output. The result would be a signal shifted by n / 2 
radians. For this model, a phase shift of -n / 2 radians was assumed. The phase shift 
can be modeled in the frequency domain as in Figure 2.5. 



a, 
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~2 



- }<*\ -> 



Figure 2.5 Phase Shift Element 



C. INPUT SIGNAL ANALYSIS 

Using the preceding development, the function of the Mach-Zender coupler can 
now be completely described. Consider an input signal, s(t), applied to an input port of 
the Mach-Zender coupler. The unused input port (2) is terminated. Let S(f) be the Fourier 
transform of s(t). The signal can be applied to a Mach-Zender coupler as shown in Figure 
2 . 6 . 
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Figure 2.6 Input Signal Through Mach-Zender Coupler 
Using Equation 2.7 the output of the first coupler becomes 
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After the delay line the inputs to the second coupler are obtained by applying the 
delay illustrated in Figure 2.4. The inputs to the second coupler are: 
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The outputs of the second coupler, b" and b [ , are obtained by applying Equation 
2.6b with 6 3 ' and b 4 as the inputs. The following result is obtained: 
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The addition of the — n 12 phase shift gives the final result in the frequency do- 
main as: 
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By taking the inverse Fourier transform of Equation 2.11 the time domain result is 
obtained: 



b"(t ) 
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s(t-T)- s(t) 
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( 2 . 12 ) 



By multiplying the matrices used to represent the phase shift, delay and two cou- 
plers (Equation 2.13), we can derive a general equation (Equation 2.14(a), 2.14(b)) where 
the output of the Mach-Zender coupler can be readily determined from the inputs S', and 
a 2 . 
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D. SUMMARY 

In this chapter the optical 2x2 coupler was shown to combine the signals of its in- 
put ports, split the signal, and send the result to the output port such that both signals 
leave the coupler simultaneously. It was determined that if only one input port was used 
then no part of the input signal was reflected back. It was also shown that the outputs 
from the Mach-Zender coupler are ( S(f)e~ j2 ” fT - S(f ))/ 2 and (S(f)e~ j2 ’* r + S(f )) / 2 
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in the frequency domain or, equivalently, ( s(t - T)- s{t)) / 2 and ( s(t -T)+ s(t )) / 2 in 
the time domain. 

In the next chapter a receiver will be proposed using the Mach-Zender coupler 
that is capable of processing a DPSK encoded signal. Using the results of Chapter II, the 
signal will be traced through the receiver and its strength relative to the receiver noise 
will be analyzed. Finally, a probability of bit-error, P [ ,, will be developed based on these 
results. 
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III. DPSK DEMODULATION ANALYSIS 



A. INPUT DPSK SIGNAL 

In Chapter I, an optical receiver capable of decoding a DPSK encoded signal was 
proposed (Figure 3.1). Chapter II then developed a mathematical model (Equation 2.12) 
of the Mach-Zender coupler to be used in this optical receiver. With these results, the 
next step is to analyze the input of a DPSK signal and assess the receiver performance 
versus the input signal-to-noise ratio (SNR). 




Figure 3.1 Optical DPSK Receiver 

A given DPSK signal will be in one of two forms (over a period of 2T where T is 
the period of the input signal) depending on the information bit being transmitted. For a 
bit “1” being sent where 0 < t < 2T , the signal will be: 

5,(0 = ±Ap T (t - T)cos(27rf c t + d 0 ) + Ap r (t)cos(2nf c t + 0 O ) . (3.1) 

For a bit “0” being sent where 0 < t < 2T , the signal will be: 

S 0 (t) = ±Ap r (t - T) cos(2nf c t + 0 O ) ± Ap T (t ) cos(2 nf c t + 6> 0 ) (3.2) 

where f c is a multiple of 1/T and p T (t — iT) with i an integer is defined as: 
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(3.3) 



jl iT < t < (i + \)T 
Pr( l iT) |q otherwise. 

Due to the symmetry of the signal, an analysis of an input 1 or 0 will produce the 
same results. To simplify the calculations an input of 0 (Equation 3.2) will be considered. 



B. RECEIVER ANALYSIS 



1. Output of Mach-Zender Coupler 

In Chapter II, the output of the Mach-Zender Coupler was determined in Equation 
2. 12. Applying Equation 2.12 to the input signal, given the transmission of a 0 (Equation 
3.2), the following output is obtained where T < / < 2T: 



« 2 (0 
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MO. 


“ 2 


± A cos(2 nf c t + 9 0 ) ± A cos(2 7rf c t + 6 0 ) 



(3.4) 



Simplifying this expression yields: 



« 2 ( 0 

MO 

which are the outputs of the Mach-Zender coupler given an optical DPSK input signal. 



± Acos(2nf c t + O 0 ) 



(3.5) 



2. Photodetector Analysis 

The next stage of the receiver is an optical photodetector (PD) which converts the 
optical signal into an electrical current signal. The PD provides an output current equiva- 
lent to the input peak power (P=A^/2) multiplied by the responsivity (93) plus the shot 
noise, r| s (t) [Ref. 4]. 

Taking our input from the Mach-Zender coupler (Equation 3.5) and applying it to 
our PD gives the following signal in the receiver over the interval where T < t < IT : 
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(3-6) 



~a 3 (t)~ 




MO 


MO. 




JR(A 2 / 2)+r] Sh (t)_ 



It is important to note at this stage that the shot noise produced by the two PD de- 
vices ( rj s (t ) and rj h (t) ) are independent of one another. 

3. Amplifier Analysis 

The next stage of the receiver consists of a transimpedance amplifier with gain Zq. 
Along with a gain to the input signal, this amplifier will introduce additive white gaussian 
noise (AWGN), r/amp(0> to the signal. In order to simplify the understanding of noise in 
the receiver a summation of total noise will be necessary. 

Both the amplifier noise, r] am p(t), and the shot noise, tj^ (t) and rj h (/) , are mod- 
eled as AWGN with zero mean. If the total noise in each individual channel is rj(t) then: 

n(0 = rj amp (t) + rj s (t) . (3.7) 

Therefore, the noise content of each channel following the amplifier stage will be rj a (t ) 
and T] h (t) with a power spectral density (PSD) of Nq. Again, the noise in each channel 
will be independent of the noise in the other channel. 

Given this simplifying assumption, the output of the amplifier is simply the input 
(Equation 3.6) plus the amplifier noise ( rj amp ( t ) ) multiplied by the amplifier gain (Zq): 



a*(t) 




Z M0 


MO. 




_Z 0 (91A 2 / 2 + rj h (t))_ 



(3.8) 



4. Integrator Analysis 

The next step in the analysis is to take the amplifier output (Equation 3.8) 
and run it through the integrator. This will give us: 
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(3.9) 
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which will simplify to: 
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(3.10) 



To continue with the analysis of the integrator output the statistical nature of the 
signal at this point will be examined. Signal a 5 (t) is defined as follows: 



fl 5 (0 = «s + N a(0 



(3-11) 



T 

where a 5 = 0 andA a (/) = J[Z 0 %(/)]<# . Since iV a (/) consists of a constant multiplied 

o 

by zero mean AWGN, its mean is also zero. Therefore a 5 represents the mean and N a (t) 
represents the noise. 

Similarly if b 5 (t ) is defined as: 



b 5 (t) = b 5 + N h ( 0 



(3.12) 



T 

where b 5 = Z 0 9L4 2 T / 2 and N h (t)= J[Z 0 ^ (t)]dt , then b 5 represents the mean and 

o 

N h (t) represents the noise. 

Next, for a gaussian random variable such as our zero mean AWGN defined by 
N a (t) and Nfj(t), the variance ( cr 2 ) is defined by: 

a N a 2 = a N„ 2 = ZqNqT. (3.13) 
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Now that the mean and variance of a 5 {t ) and b 5 (t) has been determined, an 
analysis of the decision process to choose the signal sent can be undertaken. 

5. Decision Variable Analysis 

With the mean and variance of the signals to be input to the decision element de- 
termined, an examination as the probability of an error in signal reception can begin. A 
symbol “0” was assumed to be sent to our device initially, therefore an error will have 
occurred if the symbol “1” is chosen as the output. This probability of choosing a “1” 
when a “0” is sent, Pr{ 1 10}, will be equivalent to the magnitude of the signal in channel 
a being greater than the signal in channel b given that a “0” was sent or Pr{a 5 >& 5 | 0 }. 

If a new variable, Y, is defined where: 

Y =b 5 -a 5 (3.14) 



then it is apparent that: 



Pr{a 5 > 6 5 |o} = Pr{r < 0|0} . (3.15) 

From the analysis of the integrator output it was determined that a 5 and b 5 are 
both independent gaussian random variables, and therefore Y, as the difference of a 5 and 
b 5 , is also a gaussian random variable. Consequently the mean of Y or Y can be ob- 
tained as follows: 



Y =b 5 -ci 5 =Z 0 WA 2 T/2. (3.16) 

Similarly the variance of Y(a Y 2 ) can also be obtained as follows: 

a Y 2 = a N 2 +a N 2 = Z 0 2 N 0 T + Z 2 N 0 T = 2Z 2 N 0 T . (3.17) 

Pr{l|o} or Pr{F < 0|0} can be defined by the probability density function of Y as 
follows: 
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0 

Pr{y < 0|0} = \f y (y\0)dy (3.18) 

-oo 

where the probability density function, /^(>'|0) , for a gaussian random variable is defined 
as: 



fly |o) = 



ylhti 



TIG, 



-exp 



2cr 2 



(3.19) 



Substituting Equation 3.19 into Equation 3.18 gives: 



U 0 i 

Pr{r<0|o}= \f y {y\v>)dy= J -j == — exp 

-iV2/ro- 



2o] 



dy. 



(3.20) 



(y-Y) dy 

If a new variable x is defined such that x = then dx = - — . Using this 

new variable to substitute into Equation 3.20 gives: 



YU j, 

Pr { K < 0|0} = —f== \ exp 

\2?ra v ' 



x 

T 



1 

(— G y dx) = -j== J exp 



yl2n 



YUr. 



X 

T 



dx. (3.21) 



Y 

Equation 3.21 is now in the same format as the Q function, Q(z ) , where z = — . 



Therefore: 






Pr{l|o} = Q\ — 



= Q\ 



'yj 



Z.J2NJ . 






(3.22) 
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C. PROBABILITY OF BIT ERROR 



The probability of bit error, P b , for this receiver can be determined as: 

P„ = Pr(l) x Pr{0|l} + Pr(0) x Pr{l|o} . (3.23) 

If a “1” or a “0” is equiprobable (i.e., Pr(l) = Pr(0) ) and, if by symmetry Pr{0|l} 

= Pr{l|o} , then substituting Equation 3.22 into Equation 3.23: 



P h = Q 



l , 

-9 \A 2 

V2 




(3.24) 



Further, if the peak power (P) is defined as P=A^/2 then Equation 3.24 will simplify to: 



p„=Q 



9 W 




Defining the voltage signal-to-noise ratio (SNR) as: 



Y 1 

SNR = — = -9 \A 

cr 2 



2 




(3.25) 



(3.26) 



then P h becomes: 



P„=Q(SNR). (3.27) 

Figure 3.2 shows this dependence, which illustrates the rapid decline in bit error rate as 
the signal-to-noise ratio exceeds 4 dB. 

D. SUMMARY 

In this chapter a mathematical model for a single user DPSK encoded optical re- 
ceiver was developed using the analysis of the Mach-Zender coupler from Chapter II. 
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SNR(dB) 

Figure 3.2 Pb vs. SNR 

From this model a relationship of bit-error rate versus signal-to-noise ratio was then de- 
veloped. 

In Chapter IV, a mathematical model of a multi-user CSK optical receiver will be 
developed in a very similar manner to the development in this chapter for a DPSK re- 
ceiver. Again as in this chapter, that model will be used to determine the performance 
(i.e., Pb vs. SNR) of the CSK multi-user receiver. 
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IV. CSK DEMODULATION ANALYSIS 



A. INTRODUCTION 

The preceding chapter developed a probability of bit error for a DPSK demodu- 
lated signal. A DPSK demodulated signal implies that only one user can send in informa- 
tion in the optical fiber at a time. In order to allow multiple users to utilize the fiber si- 
multaneously, the signal demodulation must be modified. As was introduced in Chapter I, 
a Code Shift Keyed (CSK) signal employing Hadamard- Walsh coding can be used to al- 
low multiple users on a fiber simultaneously. 

The Hadamard-Walsh code is a straightforward method of creating multiple or- 
thogonal signals. The Hadamard-Walsh basic 2x2 matrix is: 



Observe that the correlation of the first row vector with the second row vector is zero. 
The matrix of Equation 4.1 can be built upon to create additional orthogonal vectors as 
follows: 



1 1 

1 -1 

1 1 : 

1 -1 



-1 



1 1 
1 -1 

1 1 
1 -1 



(4.2) 



It is apparent that the individual elements of the 2x2 Hadamard-Walsh matrix are used to 
multiply the entire 2x2 matrix and Equation 4.2 can be written as: 



1111 
1 - 11-1 
11 - 1-1 
1 - 1-1 1 



(4.3) 
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Again observe that each row vector is orthogonal (i.e., has a zero correlation) to all oth- 
ers. 

The 4x4 matrix of Equation 4.3 consists of four orthogonal vectors and can be ap- 
plied to a system in which two bits or two users’ information must be decided simultane- 
ously. It is this matrix that will be used to design an optical receiver and develop the 
probability of bit error for two users on the same fiber simultaneously. 

The design of the receiver will be the same as that proposed in Chapter I and is 
shown in Figure 4. 1 . 



B. INPUT CSK HADAMARD- WALSH ENCODED SIGNAL 

A given two-user CSK Hadamard-Walsh encoded signal, S(t), will be comprised 
of the sum of two orthogonal signals, one each from channel A and channel B. Therefore: 

S{t) = S A {t) + S B {t). (4.4) 

S(t) will exist over a bit time of T—4T' where T' is the period of the input signal. Fur- 
ther: 



r = - (4.5) 

4 

where T was defined in chapter three as the period of a DPSK signal. The possible signals 
from channels A and B where 0 < t < 4 T' are: 

S ]B {t) = +Ap r ( t - 3 T') cos(2 7if c t + 0 O ) + Ap r { t - 2 T') cos(2 nf c t + 0 O ) 

+ Ap r ( t-T ') cos{27if c t +0 O ) + Ap r (0 cos(27f c t + 0 O ) (4.6) 

S ]A ( t ) = +Ap r ( t - 3 T') cos(2 nf c t + 0 O )~ Ap r ( t - 2 T') cos(2 nf c t + 0 O ) 

+ Ap r {t-T') cos(2^f c r + 0 O ) - Ap r (r) cos(2 nf c t + 0 O ) (4.7) 

S 0B ( t ) = +Ap r ( t - 3 T') cos{l7jf c t + 0 q) + Ap r ( t - 2 T') cos(2 7rf c t + 0 O ) 

- Ap r {t-T') cos(2 7rf c t + 0 O )- Ap r {t) cos(2 nf c t + 0 G ) (4.8) 
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Figure 4.1 Multiple Stage Mach-Zender Coupler Receiver 



S 0 A (/) = +Ap r ( t - 3 T') cos(2 7rf c t + 6 0 ) - Ap r (/ - 2 T') cos(2 Jif c t + 0 O ) 
- Ap r (t - T ') cos(2 7rf c t + 6 0 ) + Ap r (/) cos(2 7tf c t +8 0 ) 



(4.9) 



where the S subscript refers to the bit transmitted and the channel it was transmitted from, 
respectively. For example, Sia (0 refers to the channel A component of S(t) in which a 1 
was transmitted for channel A, and Sob (0 refers to the channel B component of S(t) for 
which a zero was transmitted for channel B. The p r (t - iT') term is defined as it was in 

Chapter ID (Equation 3.3). The frequency f c is a multiple of 1/T'. 

Again due to the symmetry of the signal an analysis of all possible inputs is not 
necessary, as each will produce the same result. Also, as each signal in channel A is or- 
thogonal to all signals in channel B, there will be no interference between channel A and 
channel B signals in S(t). To aid in the development of a probability of bit error, Pb, an 
input of a 1 in channel A and a 0 in channel B will be considered. This input corresponds 
to Equations 4.6 and 4.7 respectively and the input signal will be: 

S(t)=S 1A (t)+S 0B (t). (4.10) 



C. RECEIVER ANALYSIS 



1. Output of First Stage Mach-Zender Coupler 



In this multistage receiver the first stage is identical to that analyzed in Chapter HI 
(Figure 3.1), substituting T' for T. Again as in Chapter HI the output of the first stage is 
determined by applying Equation 2.12 to the input signal, Equation 4.9. Given a 1 in 
channel A and a 0 in channel B the following first stage output is obtained where 
T' <t < AT' : 



'SAD' 


= — cos(2 itf r t+0 n ) 


2 Ap r ( t - 3 T') + 2Ap r ( t - IT') - 2 Ap r ( t - T') 




2 


+ 2Ap r {t-3T')-2Ap r {t-T') 



A detailed derivation of this output is contained in Appendix A. 
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2. Output of Second Stage Mach-Zender Couplers 

The second-stage couplers are designed with a 2T' delay variable and conse- 
quently the analysis of their outputs is different than that of the first stage. However, with 
only a slight modification, a variation of Equation 2.12 can be used to describe their out- 
puts. Specifically for a delay of 2T' Equation 2.12 can be rewritten as follows: 



’«(*)" 


_ i 


's(t-2T')-s(t) 


m 


“ 2 


_s(t-2T') + s(t)_ 



(4.12) 



Applying Equation 4.12 to the input of the second stage, Equation 4.11, the fol- 
lowing output is obtained where 3 T' <t< 4 T': 







0 


s' u (t) 


_ i 


1 / 2(4 A cos(2 vf c t + 0 O )) 


s'oA t) 


~ 2 


l/2(-4Acos(27Tf c t + 0 o )) 






0 



which can be simplified to: 



(4.13) 



S()A (0 




0 


s; A (t) 




Acos(27if c r + 0 O )) 


S' 0B (t) 




-Acos(2^f c r + 0 O )) 






0 



These equations represents the output at the second stage Mach-Zender couplers given an 
optical multiplexed CSK Hadamard-Walsh encoded signal. 

2. Photodetector Analysis 

The remainder of the receiver for each channel is identical to the DPSK receiver 
analyzed in Chapter HI; therefore, the analysis will be similar. In the next stage of the re- 
ceiver the optical photodetector (PD) converts the optical signal into an electrical signal. 
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The PD provides an output current equivalent to the input power (P=A 2 /2) multiplied by 
the responsivity, SR , plus the shot noise, T| s (t). 

Taking our input from the Mach-Zender couplers (Equation 4.13) and applying it 
to our PD gives the following signals in the receiver over the interval 3 T' < t < AT' : 









Suit) 




3t(M J /2)+n„(0 


So, it) 




«(A 2 /2)+ij s „«) 






'U.W 



It is important to again note at this stage that the shot noise produced by the all the 
PD devices are independent of one another. 

3. Amplifier Analysis 

The next stage of the receiver, as was the case in Chapter IQ, consists of a tran- 
simpedance amplifier with gain Zq. Along with a gain to the input signal, this amplifier 
will introduce AWGN, rj amp (t), to the signal. As was the case in the analysis of the DPSK 
receiver the total noise in each channel will be summed. Since both the amplifier noise, 
f] amp(t), and the shot noise, T] s (t), are AWGN with zero mean, the total noise in each indi- 
vidual channel is T](t) and : 



ri(t) = ri amp (t) + V s (t ) . (4.16) 

The noises, tIoaU), r\iA(t), f] 0 B(t), and riiB(t), all have a power spectral density (PSD) of N 0 
and the noise in each channel is independent of the noise in the other channels. 

Given this simplifying assumption the output of the amplifier is simply the input 
(Equation 4.14) plus the amplifier noise, r\ amp (t), multiplied by the amplifier gain, Z 0 and 
can be written as: 
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(4.17) 



X'(0" 




ZoVoaIO 


Suit) 




Z 0 (KA 2 /2 + r /lA (0) 


Soil 0 




Z 0 (9iA 2 /2 + T loB (t)) 


_Su(0_ 




Z 0 tJu 10 



4. Integrator Analysis 

The next step in the analysis is to take the amplifier output (Equation 4. 1 7) and 
run it through the integrator. This will give us: 



S'o'Xit) 

Suit) 

S'o'slt) 

Suit) 



J[ Z o*7<m (')]<* 

3 r 

\[Z 0 (*A 2 U + ^mdt 

3 T 

4 T 

\lZ 0 {KA 2 12 + ^mdt 

4 r 

J[ Z o?7is(0] dt 



4 r 



3 T 



3 r 



(4.18) 



which will simplify to: 



S’oa<?) 

Suit) 

S'oslt) 

Suit) 



4 r 

\[Z 0 % A lt)]dt 

3 r 

Z 0 MA 2 T'/ 2 + J [Z 0 tj lA (t)]dt 

3 r 

Z 0 9 L 4 2 r /2 + J[Z 0 i7 0a (/)]dr 

3 r 

4 r 

j[ Z o*7is(0]<# 

3r 



(4.19) 



To greatly simplify the analysis at this point, notice the similarity of the integrator 
output defined in Equation 4.19 to that of the DPSK receiver integrator output defined in 
Equation 3.10. Specifically, the information plus noise channels in both cases (S" A (t), 
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£'"'(/) a nd b 5 (()) are identical as are the noise-only channels (S” B " (t),S 0A (t),and a 5 (t)). 
Therefore, the statistical nature of the CSK receiver integrator output can be summarized 
as: 

Soa=S, a +N oa (4.20) 

S U =S, A +N lt (4.21) 

Sob = Sqb + N ob (4.22) 

S lB =S iB +N ]B (4.23) 



where the mean of each channel can be described as S 0A and S ]B (which are both equal 
to zero) and S lj4 = S 0B = Z 0 3iA 2 T' /2. Further, the variances of each channel have identi- 
cal characteristics, and as in Chapter ID, each channel represents a gaussian random vari- 
able. Therefore the variance in each channel can be defined as: 

o N 2 =Z 2 NJ'. (4.24) 

Now that the mean and variance for each channel has been determined, an analy- 
sis of the decision process can be undertaken. 



5. Decision Variable Analysis 

With the determination of the mean and variance of the signals to be input to the 
decision element, an examination of the probability of an error in signal reception can be- 
gin. Due to the symmetry of the signal the probability of a bit error in channel A, Pb{A}> 
is equal to the probability of a bit error in channel B, Pb(B}. Therefore an error in channel 
A will be considered and the result will be the same for channel B. If a symbol 1 was 
sent in channel A, an error will occur if the symbol 0 is chosen. This probability of 
choosing a 0 when a 1 is sent in channel A, Pr^ {o|l}, will be equivalent to the magnitude 

of the signal Sqa being larger than the signal Sia . A new variable is defined as follows: 
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(4.25) 



Y = S lA -S 0A . 



It is apparent that to make an incorrect decision Y must be less than zero and therefore: 

Pr, Pr„{s„„ > Sjl}= Pr> < 0 |l}. (4.26) 



In Chapter DI a probability of bit-error was developed (Equations 3.16 to 3.22) 
relating one information plus noise channel to one noise-only channel. Similarly, in this 
CSK receiver, channel A is almost identical to the result derived in Chapter in, with the 
only difference being that the period in channel A is now T' vice T. Using the result ob- 
tained in Equation 3.24 and substituting 7' for T gives: 



Pr,W=S 



V 



— 9t<4 2 
2 




(4.27) 



Now using the definition of T', Equation 4.5, and rewriting in terms of T: 



Pr„W=<2 



-9t4 ! 



4x2N n 



(4.28) 



D. PROBABILITY OF BIT ERROR 

The probability of bit error for this receiver can be determined as: 

n=nW=nW- (4.29) 

Since the probability of an error in channel A equals the probability of an error in channel 
B, P b can be generalized by considering channel A only. In channel A: 

P b {A}= Pr„ (1) x Pr, {0|l}+ Pr„ (0) x Pr„ {l|o}. (4.30) 
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If a 1 or a 0 in channel A is equiprobable (i.e., Pr A (l)=Pr A (0)=J/2 ) and if by symmetry 
Pr A {0\l }=Pr A { I\0} then by substituting Equation 4.28 into Equation 4.30: 






4x2N 0 



T 



(4.31) 






If the power (P) is defined as P=A /2 then Equation 4.3 1 will simplify to: 




(4.32) 



and if the signal-to-noise ratio (SNR) is defined as it was in Equation 3.25, then Pb be- 
comes: 



This result implies a 3 dB degradation in performance when compared to the result for 
DPSK, Equation 3.26. Figure 4.2 shows a plot of the theoretical performance of the 
DPSK and the CSK receivers and at bit-error rates below about 3x1 O' 3 the DPSK receiver 
is degraded by 3 dB. 

E. SUMMARY 

This chapter built upon the development of the single-user DPSK receiver from 
Chapter in and has shown that the performance of a two-user CSK optical receiver is 3 
dB less than that of a single user DPSK optical receiver. The next step in the study in the 
use of the Mach-Zender coupler in an optical receiver will be to create a computer model 
of the simple DPSK receiver and obtain performance data from the model. If the model 
performance data closely resembles the theoretical data derived in Chapter IQ then vari- 
ous components of the model can be altered to conduct a sensitivity analysis of the optical 




(4.33) 
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SNR(dB) 



Figure 4.2 Pb vs. SNR (theoretical DPSK and CSK receivers) 



DPSK receiver. Chapter V will discuss the construction of the computer model and show 
the results of the sensitivity analysis. 
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V. SIMULATION AND SENSITIVITY ANALYSIS 



A. SIMULATION 



1. Design Approach 

The third chapter developed the theoretical probability of bit error for a DPSK 
optical signal demodulated using the Mach-Zender coupler. If a simulation of the optical 
and electrical receiver components can be constructed using MATLAB code and the sta- 
tistical outcomes of the simulated receiver at various signal-to-noise ratios accumulated, 
then a comparison of the simulated performance to the theoretical will confirm the accu- 
racy of the model. 



2. Mach-Zender Coupler Simulation Construction 

Previous work at the Naval Postgraduate School [Ref. 2] has extensively studied 
the physical characteristics of the Mach-Zender coupler. Given this understanding of the 
coupler, a modularized model of. the coupler was constructed in MATLAB code. Appen- 
dix B contains a copy of the code for each module. Each element of the coupler was mod- 
eled as an individual function with varying operational characteristics. 

The optical 2x2 couplers were designed with a variable input scattering matrix. 
This allows for varying amounts of coupling of the input signal(s) into the output signals. 
In addition, a delay element was constructed capable of creating any desired delay onto an 
input signal. 

The phase shift element was modeled as another 2x2 coupler as was proposed in 
the phase shift element analysis of Chapter II. A drawback of using another 2x2 coupler 
to phase shift the signal is the 3 dB attenuation of the signal by the coupler. In order to 
balance the final output power of both channels and b" of Figure 2.6) an attenuation 

element was inserted in the non-phase shifted output. As with the other modules the 
amount of attenuation was a controlled variable. 
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3. Electrical Receiver Element Simulation 



The next stage of the proposed DPSK optical signal receiver converts the de- 
modulated optical output of the Mach-Zender coupler into an electrical signal for proc- 
essing to obtain a digital output signal as described in Chapter HI. Each electrical compo- 
nent was again designed as a module with varying inputs. 

The photodetector can apply the responsivity (SR ) and add additive white gaussian 
noise to its output. The transimpendance amplifier uniformly amplifies the signal and also 
adds additive white gaussian noise. The integrator integrates the signal over one period. 
The relative strengths of the signal in channel A and channel B are then compared to 
make a decision of a “1” or a “0”. 

4. Statistical Analysis 

The MATLAB model of the receiver was repeatedly run at various SNR to obtain 
the probability of bit error for specific SNR. The simulated results and the theoretical re- 
sults are shown in Figure 5.1. The simulated results differed by about 4/10 of a dB from 
the theoretical. This indicates that the MATLAB model was accurate. 

With the accuracy of the model confirmed, the model can be used to vary the op- 
eration of various components and investigate the sensitivity of receiver performance to 
less than ideal components. Specifically, the amount of delay and the phase shift were 
altered to non-ideal values and new receiver probability of bit errors were compiled. 

B. SENSITIVITY ANALYSIS 

1. Delay Error 

In building the Mach-Zender coupler the delay element was implemented by in- 
serting an extra length of fiber equal in length to one bit period. At 100 Mbps Heinbaugh 
[Ref. 2] showed that tolerances of ±0.5% were achievable. Inserting a delay error of of 
0.5% into the MATLAB model showed no difference in Pb vs. SNR. In order to 
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Pb vs. SNR (ideal receiver) 




SNR(dB) 

Figure 5.1 Pb vs. SNR (ideal receiver) 



acquire measurable data delay, errors of 5%, 10% and 20% were introduced into the 
model. The results are shown in Figure 5.2. With 5% delay error the required SNR for a 
specific Pb increased one half of a dB. A 10% delay error resulted in over a 3 dB degra- 
dation in performance. At 20% delay error the receiver remained a random number gen- 
erator (i.e, Pb=0.5) up through a 13 dB SNR. 
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SNR(dB) 



Figure 5.2 P b vs. SNR (simulated delay error) 



1. Phase Shift Error 

The phase shift error is defined by the ability of the 2x2 coupler to change the 
phase of the input signal while keeping attenuation as close to 3 dB as possible. In the 
simulation, attenuation was maintained at 3 dB but the effects of phase delay errors of 
10° and 30° were investigated. The results are shown in Figure 5.3. 
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SNR(dB) 



Figure 5.3 Pb vs. SNR (phase shift error) 



For a 10° phase shift error receiver performance degradation was minimal and almost 
zero. For a 30° error the receiver degradation was about 2 dB. 



C. SUMMARY 

This chapter reviewed the construction of a computer model in MATLAB code to 
simulate the DPSK optical receiver. The performance of the model was compared to that 
predicted by the theoretical results obtained in Chapter in and the model performance 
closely resembled theoretical performance. Having established the accuracy of the model. 
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the model was then used to determine the sensitivity of the DPSK optical receiver given 
errors in delay and phase shift elements. Chapter VI will summarize the modeled results 
and discuss areas for future reasearch regarding the use of the Mach-Zender coupler in an 
optical receiver. 
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VI. CONCLUSION 



A. DISCUSSION OF RESULTS 

This thesis proposed the use of the Mach-Zender coupler in a single-user optical 
receiver and mathematically illustrated its theoretical implementation. Next, building on 
this implementation, the theory supporting the use of the Mach-Zender coupler to de- 
modulate multiplexed fiber-optic signals was developed and also illustrated mathemati- 
cally. From this analysis, we created a computer simulation of the single-user optical re- 
ceiver. Comparing the performance of the simulation with the expected theoretical per- 
formance established the accuracy of the model. With an accurate computer simulation 
developed, the characteristics of the model were then altered to investigate the sensitivi- 
ties of the optical receiver to non-ideal components. 

Chapter IV showed theoretically how the optical receiver could be expanded to 
demodulate a two-user multiplexed fiber optic signal. A two-user optical receiver experi- 
enced a 3-dB degradation in performance as compared to the single-user optical receiver 
proposed in Chapter IEL The receiver of Chapter IV could be extended to increase the 
number of users multiplexed on a single optic-fiber by extending the number of Mach- 
Zender coupler stages with varying delays. 

Chapter V discussed the receiver simulation using a MATLAB model, and its per- 
formance closely followed theory. This indicated that the model was accurate and that 
adjusting parameters of the model to more closely resemble an actual Mach-Zender cou- 
pler should give reasonable results in the sensitivity analysis. 

The sensitivity analysis conducted with the model showed only one-half dB deg- 
radation from a 5% delay error. As was discussed in the Chapter V, a 5% error was well 
above the 0.5% error achieved in the construction of the Mach-Zender coupler by Hein- 
baugh [Ref. 2]. Therefore, the effects of delay error due to construction on receiver per- 
formance should be minimal. 

The sensitivity analysis then varied the phase shift error up to 30° and a degrada- 
tion of less than 2 dB was experienced. A phase shift error of 30° is large; larger phase 
shift errors are unlikely. Therefore, this receiver design does not appear to be sensitive to 
probable errors in phase shift. 

This sensitivity analysis, conducted using the MATLAB model, indicates that this 
is a realizable receiver with a strong potential of simultaneously demodulating multiple 
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signals in one fiber optic line. If this method proves to be feasible, it could provide a cost- 
effective alternative to demodulating multiplexed fiber-optic signals over the current 
technique using Wavelength Division Multiplexing. 

B. FUTURE RESEARCH 

Future research on the feasibility of using the Mach-Zender coupler should focus 
on: 



1. Using the existing Mach-Zender couplers to build a DPSK optical-signal re- 
ceiver in order to confirm the sensitivity analysis already conducted with the 
MATLAB model. 

2. Constructing a CSK optical-signal receiver to investigate the reality of de- 
modulating information from two or more users simultaneously on one fiber. 
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APPENDIX A: ANALYSIS OF OUTPUT OF TWO MACH-ZENDER 



COUPLERS FOR CSK ENCODED SIGNAL 



1. SIGNAL INPUTS 

The following two signals are input to the receiver of Figure 4. 1 as S( t): 

S XA (/) = +Ap T . ( t - 3 T') cos(2 Ttf c t + 6 0 ) - Ap r ( t - 2 T') cos(2 7tf c t + 6 0 ) 

+ Ap r ( t - T') cos(2jrf c t + 6 0 )~ Ap r ( t ) cos(2 jtf c t + 6 0 ) (A. 1 ) 

S 0B (t) = +A Pr ( t - 3 r ) cos(2 7tf c t + 6 0 ) + Ap r (i f - 2 T') cos(2 nf c t + 6 0 ) 

- Ap r (t -T') cos(2 Ttf c t + d 0 ) - Ap r (/) cos(2 7tf c t + 6 0 ) (A.2) 

due to the orthogonality of the signals and superposition each signal can be 
considered separately and the end results summed. 



2. CHANNEL A INPUT ANALYSIS 



First consider the signal for channel A. To aid in analysis, cos(2 jtf c t + 0 0 ) will be 

factored out and left out of the equation until the end result is obtained. The output of 
the first coupler with a delay of T' for the channel A input will be described as S A (/) 
and S B (t). The superscript refers to the contributing input channel and the subscript 
refers to the output channel. Therefore from T' <t < AT' the output is: 

S A A (t) = ^[S IA (t-T')-S 1A (t)) = 

+Ap r (t - AT') - Ap r (t - 3 T') + A Pr (t - IT') - Ap r (t - T') 

- [+Ap r (/ - 3 T') - Ap r (t - 2 T') + Ap r ( t - T') - Ap r (r)]] 

= ^ x 2 x [-Ap r (t - 3 T') + Ap r (t - IT') - Ap r (i t-T ')] (A.3) 
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(A. 4) 



Sl{t)= X -[S 1A (t-T') + S 1A (t)} = 

\ [+Ap r (t - 47') - Ap r (/ - 37') + Ap r (/ - 27') - (/ - 7') 

+ [+4/> r (/ - 37') - zl/? r (t - 27') + (/ - 7') - Apj- (/)]] 

= — x 2 x [0] = 0 
2 

Equation A.4 is zero and can be disregarded. Therefore, applying the result of 
Equation A.3 to the input of the second stage coupler with a delay of 27' gives an output 
from 37' < t < 47' of: 

S A aA \t) = \[SA(t-T)-S A (t)} = 

^[- Ap r (/ - 57') + Ap r (t - 47') - Ap r (t - 37') 

- [-Ap r ( t - 37') + Ap r ( t - 27') - Ap r (/ - 7')]] 
= lx2x[0] = 0 (A. 5) 



SW = ^ [SA(t-T') + Sa(0] = 

i[-4p r (/ - 57') + App (t - 47') - (7 - 37') 

+ M/v <7 - 37') + zlp r <7 - 27') - zlp r (7 - 7')] 

= ^x2x [~Ap r (t - 37')] (A. 6) 

Equations A. 5 and A.6 represent the output of the second stage coupler for channel A. 
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3. CHANNEL B INPUT ANALYSIS 



Next consider the signal for channel B. To aid in analysis, cos(2 7tf c t + 6 0 ) will be 

factored out and left out of the equation until the end result is obtained. Therefore at the 
output of the first coupler with a delay of T' for the channel B input from 7' <t < 4 7' is: 



S B A (t) = UsoB(t-T')-S 0 B(t)] = 



| [+Ap r (t-4T') + Ap r (t- 37') - Ap r (t - 27') - Ap r ( t-T ') 

- [+Ap r (t - 37') + Ap r (t - 27') - A Pr (t - 7') - Ap r (f)]] 
= ^ x 2 x [-Ap r (t - 27')] 



(A.7) 



S B B (t) = ^[S IA (t-T') + S 1A (t)] = 

i [+A Pr (t - 47') + Ap r (t - 37') - Ap r (t - 27') - Ap r (t-T') 

+ [+Ap r (t - 37') + Ap r (t - 27') - Ap r (t - 7') - Ap r (f)]] 

= | x 2 x [+A Pr ( t - 37') - Ap r (t - 7')] (A.8) 



where the superscript refers to the contributing input channel and the subscript refers to 
the output channel. Equation A.7 has only one term and that term and a 27' delayed 
version of that term will be out of range of 37' < t < 47' and can be disregarded. 
Equation A.8 can be applied to the input of the second stage coupler with a delay of 2T’ 
gives an output from 37' < t < 47' of: 

S 0 B B '(t) = ^[S B (t-T')-S B (t)] = 

^[+A Pr (t - 57') - Ap r (t - 37') 

- [+Ap r (t - 37') - Ap r (t - 7')]] 

= -^ x 2 X [-Ap r (t - 37')] (A.9) 
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s!,\t)='-\s B (t-r) + s B (t)\ = 

^l+Ap r (l -ST')- Ap r (l-3T') 

+ [+Ap r ( t - 3 T') - Ap r ( t - T ')]] 

= — x 2x [0] = 0 . 

2 



(A. 10) 



Equations A.9 and A. 1 0 represent the output of the second stage coupler for channel B. 



4. OUTPUTS 

The output of the first coupler in equation 4.1 1 is as follows: 

s A (o' 

s B (0 

and consequently the output of the second coupler in equation 4.14 is as follows: 



5,"(0 + 5*(0 
5^(0 + 5j(0 



(A.l 1) 



Xo o' 




X'(o“ 


50(0 




5,1, ' ’(0 


suo 




5 0 V(0 


.5,0(0. 




.5,V(0. 



(A. 12) 



Substituting Equations A.5, A.6, A.9, A. 10 into equation A. 12 will give the result listed 
in Equation 4.14. 
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APPENDIX B: MATLAB SIMULATION PROGRAMS 



This appendix contains the MATLAB code for the various modules of the DPSK 
optical receiver model. 



1. OPTICAL DPSK RECEIVER 



% set input paramters 



f = 1e6; 

T = 1/f; 

Td = T; 

SM = [1,j;j,1]; 
L =3; 

R =1; 

No = 1; 

Zo = 1; 

PWR = [8000]; 



% frequency 
% signal period 
% length of delay 
% scattering matrix 
% attenuation (dB) 
% responsivity 
% noise PSD 
% amplifier gain 
% power levels 



t = linspace(0,T); % time intervals for one signal 



% load number of runs 
n=2000000; 



% query users for number of runs 

%n = input('enter number of simulations to run» '); 

for p=1 :length(PWR) 

A = sqrt(PWR(p)); 

% initialize data matrices 
xmtbit = zeros(1 ,n+1); 
rcvbit = zeros(1,n); 

% initial synchronization signal to receiver 

51 =A*cos(2*pi*f .*t); 
phase = 1 ; 

% create 1st data signal 
[S2,xmtbit(1)] = onechxmt(f, phase); 

52 = A*S2; 



% update phase setting 
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if xmtbit(l) == 1 
phase = -1 ‘phase; 
end 

for i = 1 :n 

% create 1 st data signal 
[S3,xmtbit(i+1 )] = onechxmt(f, phase); 

S3 = A* S3; 

% update phase 
if xmtbit(i+1 ) == 1 
phase = -rphase; 
end 

% create next input signal 
S =[S3,S2,S1]; 

% input S to Mach-Zender coupler 
[xa,xb] = mzcoup(S,T,Td,SM,L); 

% input MZ output to RCVR 
[rcvbit(i)] = rcvr(xa,xb,R,Zo,No,T); 

% progress signal to next period 

51 = S2; 

52 = S3; 

end 

% resize xmbit to remove last bit 
xmtbit = xmtbit (1 :n); 

% initialize error count 
errcnt = 0; 

% compute total nr of errors 
for i = 1 :n 

errcnt = errcnt + xor(xmtbit(i),rcvbit(i)); 
end 

Pb(p) = (errcnt/n); 

SNR(p) = 1 0*log1 0(R*PWR(p)*(sqrt(T/(2*No)))); 
end 
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% procedure to save data to .mat file 
clear P S; 
load ideal2dat; 

P=[P,Pb] 

S=[S,SNR] 

clear L A R No Zo SM t T Td i n Pb SNR errcnt SI S2 S3 rcvbit xmtbit; 
clear f PWR phase p xa xb; 
save ideal2dat; 
clear 

% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 



2. ONE CHANNEL TRANSMITTER 

function [S, xmtbit] = onechxmt(f, phase) 

% 1 CHANNEL DPSK TRANSMITTER one bit output 
% [S, XMTBIT] = ONECHXMT(F) 

% Creates a random transmit (XMTBIT) bit of one or zero. 
% Based on transmit bit it returns a signal 
% (S) one period in length that is inphase 

% if XMTBIT is a zero and 180 degrees out of 
% phase if xmtbit is one. 

% decide which bit transmitter will simulate 
% signal to send 

T=1/f; % compute period of signal 

t=linspace(0,T); % generate 100 time intervals 
x= rand - 0.5; 



if x <= 0 
xmtbit=0; 

S=phase*cos(2*pi*f.*t); 

% compute S for XMTBIT = 0 

else 

xmtbit=1 ; 

S=-1 *phase*cos(2*pi*f.*t); 

% compute S for XMTBIT = 1 
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end 



% Naval Postgraduate School, ECE 
% Created by: LT Paul W. Wehner 



3. MACH-ZENDER COUPLER 



function [a,b] = mzcoupsm(S,T,Td,SI\/l1 ,SM2,SM3,L) 

% MACH ZENDER COUPLER with varied scattering matrices output 
% For a given input signal, S (with a period T, 

% three periods in length), returns two signals, 

% A and B (both one period in length) where 
% A = S(t-T_actual)+S(t) and B=S(t-T_actual)-S(t). 

% T_ACTUAL indicates the actual length 
% (i.e., seconds) of the coupler delay 

% element and SM# indicates the 
% scattering matrix of the individual #'rd 
% 2x2 optical couplers. L is the attenuation 

% in channel a to equalizer signal strength 
% lost in channel b due to the last phase 
% shift. 

% determine # discrete quantities in 
% a signal period (T) of S 
n = length (S)/3; 

% demux signal S in 1st 2x2 coupler 
[al ,b1 ] = optcoup2(SM1 ,S); 

% delay signal in channel a 
[a2,d] = delay(a1,T,Td); 

% adjust length of channel a and 
% channel b signals 
a3 = a2(d-n+1:d); 
b2 = b1(d-n+1:d); 

% 2nd 2x2 coupler output 
[a4,b3] = optcoup2(SM2,a3,b2); 

% phase shft element output 
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[q,b4] = optcoup2(SM3,b3); 

% attenuate channel a 
[a5] = atten(a4,L); 

% define MZ coupler output 
a = real(a5); 
b = real(b4); 

% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 



4. OPTICAL COUPLER 



function [a,b] = optcoup2(SM,S1,S2) 

% optical 2x2 coupler 
% OPTCOUP2(SM,S1,S2) 

% For a given input signal, SI , and a possible 
% second signal S2, provides two outputs, 

% A and B. Both outputs are attenuated by a 
% factor of 3dB and multiplied by the scattering 
% matrix, SM. SM must be [2x2]. 

% determime if second signal S2 was input 
if nargin == 2 
S2=zeros(1,length(S1)); 
elseif length(SI) ~= length(S2) 
disp('ERROR - Input arrays (SI and S2) must 1 ) 
disp(' be of the same length') 
a=0; 
b=0; 
return; 
end 

% define input matrix 
y = [S1;S2]; 

% multiply by scattering matrix 
x = SM*y; 

% attenuate output by 3dB 
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x = x./sqrt(2); 



% define individual outputs 
a = x(1 
b = x(2,:); 



% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 



5. DELAY ELEMENT 



function [a,d] = delay (a,T,Td) 

% DELAY ELEMENT 

% DELAY(A,T,Td), for a given signal A (A must 
% be three periods in length with a period of 
% T seconds) outputs a new signal A delayed Td 
% seconds. Indicates at which element 
% value, N, in a the delayed vector begins 
% where the vector is listed right to left 

% create working copy of a 
b = a; 

% determine # discrete quantities of one 
% period of a 
i = length(a)/3; 

% compute delay as percentage of actual delay 
perc_Td = Td/T; 

% equate percent delay to # of discrete 
% quantities of a 

j = round(i * perc_Td); % quantized delay 
k = (3*i) - j; % new start of a 



% zero out delayed discrete elements of a 
for I = (3*i):-1:(k+1) 
a(l)=0; 
end 
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m = 3*i; 



% overwrite a with delayed version of itself 
for n= k:-1:(2*i-j+1) 
a(n)=b(m); 
m=m-1; 
end 

d=k; 

% Naval Postgraduate School, ECE 
% Created by: LT Paul W. Wehner 



6. ATTENUATOR 

function [a] = atten(a,L) 

% ATTENUATOR 
% ATTEN(A,L) 

% For given input signal, A, performs 
% attenuation on signal equal to 
% L dB. 

a = a / (sqrt(1 0 A (L/1 0))); 

% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 



7. ELECTRICAL CIRCUIT ELEMENTS 

function [rcvbit] = rcvr(xa,xb,R,Zo,No,T) 

% OPTICAL DPSK SIGNAL RECEIVER OUTPUT 
% [RCVBIT]=CRCVR(XA,XB,R,ZO,NO,T) 

% For given input signals, XA and XB 
% this applies photodetection 
% proportional to the responsivity (R) to 
% convert the optical signal to electrical, 

% amplifies the electrical signal 
% with a gain of (Zo) and then 
% integrates the signal over a 
% single time period (T) to obtain signal 
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% energy in each branch. A comparison of 
% two branches is then made to 
% determine the transmitted bit. Noise is 
% added in the photodetector and 
% the amplifier with a power spectral 
% density of No/2. 

% inut optical signal to photodetector 
xal = pd(xa,R,No,T); 
xbl = pd(xb,R,No,T); 

% input signal to transimpedance amplifier 
xa2 = tranamp(xa1 ,Zo,No,T); 
xb2 = tranamp(xb1 ,Zo,No,T); 

% integrate both channels over T 
xa3 = intg(xa2,T); 
xb3 = intg(xb2,T); 

% compare channels to decide received bit 
if xa3 > xb3 
rcvbit = 1 ; 
else 

rcvbit = 0; 
end 

% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 



8. PHOTODETECTOR 

function [a] = pd(x,R,No,T) 

% PHOTODETECTOR 
% PD(x,R,No,T) 

% For given input signal, X, performs 
% photdetection on signal power proportional 
% to the responsivity, R, in amps/watt. 

% Assume signal is output over a 1 ohm 
% load such that amps correlate directly to 
% volts. Adds photodetection shot noise to 
% output signal A with PSD(No/2) as defined 
% from user input. 
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% Create shot noise 
s = randn(1 ,length(x)); 

% adjust shot noise PSD 
amp = sqrt(No/T); 
s = amp * s; 

% compute signal power 
Px = (x. A 2)./2; 

% multiply input power by responsivity 
% to get output in amps 
Ax = R * Px; 

% add shot noise to create output signal 
a = Ax + s; 

% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 



9. TRANSIMPEDANCE AMPLIFIER 

function [a] = tranamp(x,Zo,No,T) 

% TRANSIMPEDANCE AMPLIFIER 
% TRANAMP(x,Zo,No,T) 

% For given input signal, X, performs 
% uniform amplification on signal 
% with gain Zo. Adds noise to signal prior 
% to amplification with PSD No/2. 

% Assume signal is output over a 1 ohm 
% load such that amps correlate directly to 
% volts. Adds photodetection shot noise to 
% output signal A with PSD(No) as defined 
% from user input. 

% Create amp noise 
s = randn(1 ,length(x)); 



% adjust amp noise 
amp = sqrt(No/T); 
s = amp * s; 
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% add amp noise to input signal 
a = x + s; 

% adjust output by gain, Zo 
a = Zo * a; 



% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 



10. INTEGRATOR 

function [a] = intg(x,T) 

% INTEGRATOR 
% INTG(x,T) 

% For given input signal, X, performs 
% trapazoidal integration on signal 
% over one time period T. Assumes 
% elements of X are equally spaced 
% from 0 to T. 

a = (T/length(x)) *trapz(x); 



% Naval Postgraduate School 
% Created by: LT Paul W. Wehner 



11. Q FUNCTION 
function [Q,P] = qfn(x) 

% This function computes the Gaussian Q-function 
% using the rational approximation 26.2.17 of 
% Abromowitz and Stegun. For sufficiently large 
% arguments, the asymptotic expansion 26.2.12 is used. 
% The rational approximation is accurate to within 
% 7.5e-8. Q is the value of the Q-function (integral of 

% Gaussian pdf) of zero mean and unit variance from the 
% input quantity to infinity, and P = 1 - Q is the 
% cumlative distribution function. The call is 
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% [Q,P] = qfn(x). 

% 

% R.E. Ziemer 
% 4/17/94 

b(l)=0.31938153; 
b(2)=-0. 356563782; 
b(3)=l. 78 1477937; 
b(4)=- 1.82 1255978; 
b(5)=l. 330274429; 
p=0.2316419; 
y=abs(x); 

T=l/(l+p*y); 

TT=zeros(l,5); 
for r= 1:5 
TT(r)=T A r; 
end 

Z=exp(-y A 2/2)/sqrt(2*pi); 
if y<4 

Ql=Z*b*TT; 

else 

Ql=(Z/y)*(l-l/y A 2+3/y A 4); 
end 
if x>0 

Q=Qi; 

else 

Q= l -Q l ; 

end 

P=i-Q; 
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